Collapsing transition of spherical tethered surfaces with many holes 



Hiroshi KoibuchQ 

Department of Mechanical and Systems Engineering, 
Ibaraki National College of Technology, Nakane 866 Hitachinaka, Ibaraki 312-8508, Japan 

(Dated: February 6, 2008) 

We investigate a tethered (i.e. fixed connectivity) surface model on spherical surfaces with many 
holes by using the canonical Monte Carlo simulations. Our result in this paper reveals that the model 
has only a collapsing transition at finite bending rigidity, where no surface fluctuation transition can 
be seen. The first-order collapsing transition separates the smooth phase from the collapsed phase. 
Both smooth and collapsed phases are characterized by Hausdorff dimension H ~ 2, consequently, 
the surface becomes smooth in both phases. The difference between these two phases can be seen 
only in the size of surface. This is consistent with the fact that we can see no surface fluctuation 
transition at the collapsing transition point. These two types of transitions are well known to occur 
at the same transition point in the conventional surface models defined on the fixed connectivity 
surfaces without holes. 

PACS numbers: 64.60.-i, 68.60.-p, 87.16.Dg 



I. INTRODUCTION 

Two-dimensional surfaces emerge as interfaces between 
two different materials. Biological membranes are consid- 
ered as such two-dimensional interfaces separating some 
biological materials. The model of Helfrich, Kleinert, and 
Polyakov (HPK) [EH, 3 describes the mechanics of such 
surfaces on the basis of the two-dimensional differential 
geometry and the statistical mechanics [1, [H, H| . 

One interesting topic is the crumpling transition, 
which has long been inve stig ated theoretically B, H, Q 
and numerically [H El El El El EE El IlJ The 
surface fluctuations grow anomalously at certain finite 
bending rigidity in the surface models, where a surface 
collapsing phenomena can also be seen at the same transi- 
tion point. These two phenomena have been understood 
as the crumpling transition. The surface fluctuation of 
the model is characterized by fluctuations of the bend- 
ing energy S2 and, consequently, we can see an anoma- 
lous peak in the corresponding specific heat Cs 2 ■ Cm the 
other hand, the collapsing phenomena are characterized 
by a discontinuous change of the mean square size X 2 
and, then, we can see a discontinuous (or a continuous) 
change in the roughness exponent or in the Hausdorff 
dimension. 

We have not yet seen that any surface models undergo 
only one of the two transitions. Experimentally, it was 
reported that a surface fluctuation transition is accompa- 
nied by a collapsing transition in certain artificial mem- 
branes [l9[ . However, it is possible to consider that these 
phenomena are two different ones. Buckling of a thin 
elastic shell is one of the collapsing phenomena [2(| , and 
it is not always accompanied by the surface fluctuations. 
In biological membranes, the surface fluctuation called a 
rippling transition [2l| is not always accompanied by the 
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surface collapsing phenomena. 

In this paper, we show that the collapsing transition 
occurs in a fixed connectivity surface model at finite 
bending rigidity b c , where no surface fluctuation tran- 
sition is observed. The surface model is defined by the 
standard HPK Hamiltonian, which is discretized on tri- 
angulated spherical surfaces with many holes. The start- 
ing surface configuration for the Monte Carlo (MC) sim- 
ulations is constructed such that the ratio of the area of 
the holes to that of the surface including the holes is fixed 
in the thermodynamic limit. 

The transition point b c (— 1.4 ~ 1.5) is relatively larger 
than that (6 C ~0.77) in the same model on spherical sur- 
faces without the holes. As a consequence, the surfaces 
with holes are relatively smooth at the collapsing transi- 
tion point b c . Then, we have H~2 even in the collapsed 
phase. 

The result in this paper indicate a possible collapsing 
transition which is not accompanied by the surface fluc- 
tuation transition in biological or artificial membranes, 
although the self-avoiding property [H HI H3| is not as- 
sumed in the model. A crumpling phenomenon can also 
be seen on thin elastic sheets and has been investigated 
by singularity analysis [H[26j]. However, the collapsing 
transition in biological membranes is not yet completely 
understood because of the complexity of biological mem- 
branes. 



II. MODEL AND MONTE CARLO TECHNIQUE 

In order to construct the lattice that has holes, we 
start with the icosahedron. Firstly, we divide the edges 
of the icosahedron into £ pieces of the uniform length a. 
Then, we have a triangulated surface of size No = 10£ 2 +2, 
which is the total number of vertices in the triangu- 
lated sphere without holes. Then an edge length of the 
icosahedron corresponds to I edges in the triangulated 
sphere. Secondly, the £ edges are divided into m groups 
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(to = 1, 2, • • ■ ), where each group has length L — l/m in 
the unit of a if to divides t. As a consequence, we have 
a sublattice of edge length La in the triangulated sphere 
of size Nq. We should note that this sublattice makes 
compartments on the lattice. Finally, one part of the 
compartments in the sublattice is labeled as holes, and 
the other remaining part is labeled as the lattice points. 
Thus, we have a triangulated lattice with many holes. 
The total number of holes in one face of the icosahedron 
is given by m(m — l)/2 and, hence the total number of 
holes is 10m(TO— 1) on the lattice. Because of the holes on 
the surface, the total number of vertices N of the lattice 
with holes are reduced from No, which is the total num- 
ber of the original triangulated lattice as stated above. 
We note that N includes the vertices on the boundary of 
holes. The size of lattice can be characterized by (A, to). 




(a) (N,m) = (1892, 3) 



(b) (N,m) = (2202, 4) 



FIG. 2: (Color online) Starting configuration of surfaces of 
(a) (N,m) = (1892, 3) and (b) (N,m) = (2202, 4). The total 
number of holes is 60 in (a), and it is 120 in (b). 




(a) m = 3 



(b) m = 4 



FIG. 1: A face of icosahedron corresponding to (a) m — 3 
and (b) m = 4. The shaded triangles excluding their corners 
correspond to the holes. The total number of holes in the face 
is (a) 3 and (b) 6, which are given by m(m — l)/2. 

We show the holes in a face of the icosahedron in 
FigsOJa) an d Hlb), which respectively correspond to 
to — 3 and m = 4. Note that the faces are those of the 
triangulated surfaces, whose total number of vertices is 
given by the above described number Ao = 10£ 2 +2. The 
small triangles including the shaded ones in the figures 
form the above-mentioned sublattice or compartments 
and, therefore, they are triangulated with more fine tri- 
angles. At the corners of each shaded triangle, the three 
fine triangles, which are not shown in the figures, are 
excluded from the hole (or equivalently included in the 
surface) . 

Figures Ufa) and EJb) are the lattices of (A, to) = 
(1892,3), and (N,m) = (2022,4). We note that the ra- 
tio R m (N) of the area of holes to that of the surface 
including the holes is R m (N)—>-(m — l)/2m in the limit 
of N — > oo or equivalently Ao — > oo. In fact, the total 
number of triangles in a hole is L 2 — 3 and, then, the to- 
tal number of triangles in the holes is 10to(to— 1)(L — 3), 
which is easily understood since the total number of faces 
in the icosahedron is 20, and the total number of holes in 
a face is m(m— 1)/2 as stated above. On the other hand, 
the total number of triangles on the triangulated sphere is 
27V -4. Then, we have R m = 10to(to-1)(L 2 -3)/(2A -4). 



By using L 2 = (t/m) 2 and i 2 = (N a 
Rm = {tti — 1)/2to in the limit of Nq — 



- 2)/ 10, we have 
oo. 



We assume two values for to such that 



to = 3, to = 4, 



(1) 



and four sizes of lattices are assumed for each to. Table 
U shows the total number of vertices N and Ao for each 
to and the corresponding R m (N), where Ao is the total 
number of vertices of the original lattice, which has no 
holes. We find that the ratio R m (N) in Table U are very 
close to the values #3(00) = 1/3 or #4(00) = 3/8, although 
A is finite in Table |H 



TABLE I: The total number of vertices N and No for each 
m, and the corresponding R m (N). Nq is the total number of 
vertices of the original lattice, which has no holes. 

m N(N ) R 3 (N) m N(N ) R 4 (N) 



3 1892(2252) 0.3320 

3 3512(4412) 0.3327 

3 5612(7292) 0.3329 

3 8192(10892) 0.3330 



4 2202(2562) 0.3735 

4 4402(3282) 0.3736 

4 6042(7842) 0.3744 

4 7722(10242) 0.3746 



The model is defined by the partition function 

,/ N 

Z = J Y[dX lC xp[-S{X)}, 



(2) 



i=l 



S(X) = S 1 + bS 2 , 



where b is the bending rigidity, J denotes that the center 
of the surface is fixed in the integration. S(X) denotes 
that the Hamiltonian S depends on the position variables 
X of the vertices. The Gaussian bond potential Si and 
the bending energy S2 are defined by 



(ij) 



^(l-n^nA, (3) 

m 



where Yl(ij) m ^1 i s the sum over bonds (ij) connecting 
the vertices i and j including those on the boundary, and 
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Yl(ij) m ^2 is also the sum over bonds (ij), which are 
edges of the triangles i and j. in Eq. ([3]) is the unit 
normal vector of the triangle i. We note that S2 is a 
quantity that is defined on the bonds. Note also that S2 
is not defined on the boundary bonds, because rij is not 
defined inside the holes. 

The unit of b is fcT, where k is the Boltzmann constant 
and T is the temperature. Let A be the string tension 
coefficient such that S — X Si+bS 2, then the length unit a 
of the model is given by y/kT/X. Consequently, the unit 
of Si can be expressed by a 2 or kT/X. Note also that we 
can arbitrarily choose the length unit a because of the 
scale invariant property of Z in Eq.([2]). Therefore, A=l 
was assumed in the simulations. 

The total number of Monte Carlo sweeps (MCS) af- 
ter the thermalization MCS are about 5 ~ 6 x fO 8 for 
(N,m) = (8192,3), 3 - 4 x 10 8 for (JV, m) = (5612,3), 
and 2 - 3 x 10 8 for (JV, to) < (3512,3) at b close to the 
transition point. Relatively smaller number of MCS was 
performed at non-transition points. In the case m = 4, 
the total number of MCS after the thermalization MCS 
are 5 - 6 x 10 8 for (JV, to) = (7722,4), 3 - 4 x 10 s for 
(JV, to) = (6042, 4), and 2 - 3 x 10 8 for (JV, to) < (3282, 4). 

We use the canonical MC technique to simulate the 
integrations in the partition function. The vertices X are 
shifted so that X' = X+SX, where 6X is randomly chosen 
in a small sphere. The new position X' is accepted with 
the probability Min[l, exp(— AS)], where AS = S(new) — 
S (old). 

The radius of the small sphere for 8X is chosen so that 
the rate of acceptance for X is about 50%. We introduce 
the lower bound 1 x 10 -8 for the area of triangles. No 
lower bound is imposed on the bond length. 



III. RESULTS 

Firstly, we show in FigOHa) a snapshots of the (JV, m) = 
(8192, 3) surface, which is obtained in the smooth phase 
at the transition point b = 1.39. Figure ^[b) shows a 
snapshot in the collapsed phase of the same surface at the 
same transition point. The surface sections of Figs[3ja) 
and[3fb) are shown in Figs|3fc) and[3]Jd), respectively. 
We find from the surface sections that there is a lot of 
empty space inside the surface even in the collapsed sur- 
face. This implies that the surface is not completely col- 
lapsed in the collapsed phase. The vertices are not al- 
ways confined in a small region; they are considered to 
be spread over some two-dimensional region in R 3 such 
that the surface becomes almost smooth. 

The mean square size X 2 is defined by 

* 2 = ^E(^) 2 > * = ^£*. ( 4 ) 

i i 

where X is the center of the surface. The mean square 
size of the snapshots in FigJ3Ja) is X 2 = 127, whereas 
X 2 = 69 in FigOJb). These values of X 2 are typical 



(a) 6=1.39 (smooth) (b) 6 = 1.39 (collapsed) 




(c) The surface section (d) The surface section 



FIG. 3: (Color online) Snapshots of surfaces of (JV, m) = 
(8192,3) at b = 1.39; (a) a smooth surface, (b) a collapsed 
surface, (c) the surface section of (a), and (d) the surface 
section of (b). The mean square size X 2 = 127 in (a), and 
X 2 = 69 in (b), where X 2 is defined by Eq.©. figures were 
drawn in the same scale. 




Oi 1 L 1 1 1 1 — 1 

1.3 1.35, 1.4 1.4 1.5 

(a) b (b) b 

FIG. 4: The mean square size X 2 vs. b obtained on the 
surfaces of (a) m = 3, (b) m = 4, and (c) m = 5. 



to the smooth phase and to the collapsed phase on the 
(JV, to) = (8192, 3) surface. 

Figures Ht^a) and HJh) show X 2 versus b obtained on 
the surfaces of to = 3 and to = 4, respectively. The vari- 
ation of X 2 versus b seems to be not so sharp c omp ared 
with the results of the tethered surface model in [12| even 
when JV becomes large. The surface size seems to change 
continuously in the model in this paper. 
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FIG. 5: The fluctuation C x i vs. b on the surfaces of (a) m — 3 
and (b) m — 4. The error bars denote the statistical errors, 
which were obtained by the standard binning analysis. 



log-log scale. The straight lines drawn in the figures were 
obtained by fitting the data to 

<™ ~ N a , (6) 

where a is a critical exponent of the transition. Then, 
we have 

cr = 1.18 ±0.07 (m = 3), 

<T = 1.13 ±0.08 (to = 4). (7) 

These results indicate that the collapsing transition is of 
first order. In fact, we understand from the finite-size 
scaling theory [U [28| that the maximum C™2 X should 
scale according to N a (a = 1) at the transition point if 
the transition is of first order 1231. 



The fluctuation Cx 2 of X 2 can be defined by 

C x ^^{(X 2 -{X 2 )f), (5) 

which is expected to reflect surface collapsing phenom- 
ena. Figures[5fa) and[S](b) are C x z against b obtained on 
the lattices of to = 3 and to = 4, respectively. Sharp peaks 
can be seen in Cx 2 and grow larger and larger when the 
size N increases. This anomalous behavior of Cx? indi- 
cates a discontinuous change of the surface size. There- 
fore, this surface collapsing phenomenon can be viewed 
as a phase transition contrary to the expectation from 
the smooth variation of X 2 in Figs[Ua) andHJb). 

From Figs 0^a) and E^b) we find that the transition 
points b c are about b c = 1.4 and b c = 1.5 in the cases 
to = 3 and m = 4, respectively. Both of them are rela- 
tively larger than the transition point b c ~ 0.77 of the 
same model on the fixed connectivity surface without 
holes P3. 




2000 5000 _ 2000 5000 



FIG. 6: Log-log plots of the peak values C™2 X against N, 
which were obtained on the surfaces of (a) m — 3 and (b) 
m — A. The straight lines drawn in the figures were obtained 
by fitting the data to Eq.©. 

In order to see the order of the transition, we plot the 
peak values C^?2 X against N in Figs (Ha) and [f^b) in a 




d) x (e) '" x ^ 



FIG. 7: The variation of X 2 against MCS obtained on the 
surfaces of (a) (N,m) = (3512,3) at b = 1.34, (b) (N,m) = 
(5612,3) at 6=1.4, and (c) (N, m) — (8192, 3) at b=1.39. The 
corresponding distribution h(X 2 ) of X 2 are shown in (d), (e) 
and (f). 

Figures a), EJb), and [Tic) show the variation of X 2 
versus MCS obtained at 6= 1.34 on the (N, m) = (3512, 3) 
surface, b = 1.4 on the (N,m) = (5612,3) surface, and 
b = 1.39 on the (N,m) = (8192,3) surface, respectively. 
Figures [7jd), [7te), and [714) are the normalized distri- 
bution h(X 2 ) (or histogram) of X 2 corresponding to the 
variations of X 2 in Figs[7jja),[7lb), andlljc), respectively. 
We find a double peak structure in h(X 2 ); one of the 
peaks represents the swollen and smooth state, which 
corresponds the snapshot in FigJ3Ja), and the other rep- 
resents the collapsed state, which corresponds the snap- 
shot in FigO^b). 

Figures [5]J a), [5Jb), and [51c) show the variation of X 2 
versus MCS obtained at 6= 1.44 on the (N, m) = (3282, 4) 
surface, b — 1.5 on the (N,m) = (6042,4) surface, and 
b= 1.51 on the (N,m) = (7722,4) surface, respectively. 
Figures Eld) , [SJe) , and[5]Jf) are the normalized histogram 
h(X 2 ) of X 2 corresponding to the variations of X 2 in 
Figs[5Ja), [Sfb), and[5Ic), respectively. The double peak 
structure can also be found in h(X 2 ) in the case to = 4. 

The double peaks in h(X 2 ) allow us to estimate the 
Hausdorff dimension H in the smooth state and that in 
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FIG. 8: The variation of X 2 against MCS obtained on the 
surfaces of (a) (N,m) = (3282,4) at b = 1.44, (b) (N,m) = 
(6042, 4) at 6= 1.5, and (c) (N, m) = (7722, 4) at 6= 1.51. The 
corresponding distribution h(X 2 ) of X 2 are shown in (d), (e) 
and (f). 



FIG. 9: Log-log plots of X 2 (col) and X 2 (smo) against N, 
which were obtained on the surfaces of (a) m = 3 and (b) 
m = 4. The straight lines drawn in the figures were obtained 
by fitting the data to Eq.j8|. 



the collapsed state, where H is defined by 

X 2 oc N 2 ' H . (8) 

Note that we understand from the value of H how smooth 
the surface is, because H is defined so that H varies 
depending on the distribution of vertices. In fact, we 
have H — > 2 when the vertices form a two-dimensional 
region in R 3 , and we have H—>3 when the vertices form 
a three-dimensional region of constant density of vertices. 
Moreover, the vertices can form a zero-dimensional region 
(= a point) if no self-avoiding property is assumed, and 
then we have H — > oo in this case. 

TABLE II: The lower bound X^f° l and the upper bound 
X^ax for obtaining the mean value X 2 (col) in the collapsed 
state, and the lower bound X^f™° and the upper bound 
X^™° for obtaining the mean value X 2 (smo) in the smooth 
state. 



m 


N 


y2 col 


y2 col 


y'2 smo 


Y'2 smo 


3 


1892 


11 


21 


23 


34 


3 


3512 


21 


39 


42 


65 


3 


5612 


40 


70 


72 


108 


3 


8192 


50 


100 


110 


160 


4 


2202 


12 


24 


26 


39 


4 


3282 


22 


39 


42 


60 


4 


6042 


35 


70 


75 


118 


4 


7722 


50 


85 


88 


145 



In order to obtain X 2 in the collapsed state, we evalu- 
ate the mean value X 2 (co\) by averaging X 2 in the range 
X 2 f°^ <X 2 < X 2c °\ where X 2 f° l and X 2c ° l are shown 

Illlll IllcLX. ' Illlll II leLX 

in Tabic [TTJ The mean value X 2 (smo) in the smooth 
state can also be evaluated by averaging X 2 in the range 
X 2 jr<X 2 <XHr, where X^° and Xl»™ are also 
shown in Table HT1 

The mean values X 2 (col) and X 2 (smo) obtained on 
the surfaces m = 3 and m = 4 are respectively plotted in 



Figsjnja) and[^b). The error bars denote the standard 
deviations, which were obtained by splitting the range of 
X 2 shown in Table ITU into 10 ranges; as a consequence, 
the errors become large (small) when the range is wide 
(narrow). The straight lines in the figures were drawn 
by fitting the data to Eq.© with the weight of inverse 
errors. Then, we have Hausdorff dimensions H cal and 
jjsmo j n collapsed state and the smooth state such 
that 

H co1 = 1.86 ± 0.19, H smo = 1.87 ± 0.12 (to = 3), 
H co1 = 1.97 ± 0.17, H smo = 2.02 ± 0.25 (to = 4). (9) 

The values of H in Eq.© are all close to H = 2, which 
is the topological dimension of the two-dimensional sur- 
faces. Thus, we understand that not only the smooth 
state but also the collapsed state can be viewed as a 
smooth surface. Only difference between the smooth 
state and the collapsed state is in the size, which is char- 
acterized by X 2 . No difference can be seen in the Haus- 
dorff dimensions of the smooth state and that of the col- 
lapsed state. 



S 2 /N B S 2 /N B 




1.3 1.35, 1.4 1.4 , 1.5 

(a) D (b) b 

FIG. 10: The bending energy S2/NB versus 6 obtained on the 
surfaces of (a) m — 3 and (b) m — 4. Nb is the total number 
of bonds. 
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Figures fTOTal and [TOTb) show the bending energy 
S%/Nb versus b obtained on the surfaces of m = 3 and 
m = 4, respectively. The values of S2/NB smoothly varies 
against b close to the transition point even when N is 
increased. This indicates that there is no surface fluctu- 
ation transition. 




1.3 1.35 1.4 1.4 1.5 

(a) b (b) b 

FIG. 11: The specific heat Cs 2 vs - b on the surfaces of (a) 
m = 3 and (b) m — 4. The error bars are the statistical errors, 
which were obtained by the binning analysis. 

In order to confirm that no surface fluctuation transi- 
tion occurs, we plot the specific heat defined by 

C S2 = ^((S 2 -(S 2 )f) (10) 

in Figslllfa) and lllf b). which respectively correspond 
to S 2 /N B in FigsEEHa) and EEDfb). We find a peak in 
each Cs 2 ; and this implies that the surface fluctuation 
becomes large at the peak point. However, the peak val- 
ues C!? ax decrease with increasing N. This obviously 
indicates that the surface fluctuation phenomenon is sup- 
pressed at large N and, hence, the surface fluctuation is 
not considered to be a phase transition. 

IV. SUMMARY AND CONCLUSION 

A tethered surface model has been investigated on 
triangulated spherical surfaces with many holes by the 



canonical Monte Carlo simulations. The ratio of the 
area of holes to that of the surface including the holes 
is fixed in the starting configurations, and it is given by 
(m — 1) I Ira in the limit of N — > 00, where m(= 3, 4) rep- 
resents the number of partitions of an edge of the icosa- 
hedron. As a consequence, the total number of holes 
remains unchanged when the size N increases under the 
condition that m is fixed. 

We found a first-order transition of surface collaps- 
ing phenomena, which are characterized by a discontinu- 
ous change of X 2 . This transition separates the smooth 
phase from the collapsed phase. The surfaces are smooth 
in both phases and, then, the Hausdorff dimension in the 
smooth state and that in the collapsed state are almost 
identical, and they are H ~ 2, which is identical to the 
topological dimension of surfaces. For this reason, it is 
possible to consider that the transition is observed in bi- 
ological or artificial membranes, which are not always 
closed. The model surface can be considered as an open 
surface, because the size La of holes is comparable (i.e. 
not negligible compared) to its surface size in the limit 
of N->oo. 

Moreover, the collapsing transition is not accompanied 
by the surface fluctuation transition, which is the one 
characterized by a discontinuous change of the bending 
energy S 2 . A surface fluctuation phenomenon occurs in 
the model and, therefore, the specific heat Cs 2 has a peak 
at the surface collapsing transition point. However, the 
peak disappears from Cs 2 at sufficiently large N. 
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